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ABSTRACT 



There is now convincing evidence that the intensity of HCN molecular line emission is enhanced around active galactic nuclei. In this 
paper we examine the specific case of the Seyfert galaxy NGC 3227, for which there are subarcsecond resolution data for the HCN 
(1-0) 88 GHz and CO (2-1) 230 GHz rotational lines, enabling us to spatially separate a circumnuclear ring at a radius of 140 pc and an 
inner nuclear region within 40 pc of the AGN. The HCN(1-0)/CO(2-1) flux ratio differs by more than an order of magnitude between 
these two regions. We carry out large velocity gradient (LVG) computations to determine the range of parameters (gas temperature and 
density, HCN/CO abundance ratio, column densities and velocity gradients) that yield physically plausible solutions for the observed 
flux ratio in the central 100 pc. The observed HCN/CO intensity ratio in the nucleus is consistent with very optically thick thermalized 
emission in very dense (> 10 5 cirT 3 ) gas, in which case the HCN/CO abundance ratio there is unconstrained. Alternatively, the 
HCN/CO intensity ratio could be due to optically thinner emissions but with very high (~ 10~ 2 ) HCN/CO abundance ratios. This 
possibility is more consistent with the CO and HCN emissions seen in the nuclei of the Seyfert galaxies NGC 1068 and NGC 6951. 
It would imply the velocity gradients are large and the clouds may be gravitationally unbound. We estimate that the X-ray ionisation 
rate at radii less than 20 pc in the centre of NGC 3227 exceeds 10~ 13 s -1 . X-ray ionisation and heating may lead to high HCN/CO 
ratios in warm gas in a high-ionisation molecular phase near the AGN. 

Key words. Galaxies: active - Galaxies: individual (NGC3227) - Galaxies: ISM - Galaxies: nuclei - Radio lines: galaxies 



1. Introduction 

During the last decade there has been an increasing observa- 
tional effort to understand the nature of the HCN emission from 
galaxies. This i ncludes studi e s of HCN emi s sion in nearby AGN 
and s t arbursts dKohno et all I200U iKohnol I2005I: iKohno et all 
l2008t iKrips et all I20071 l2008l) as well as luminous and ultra- 
luminous galaxies ( Gracia-Carpio et all l2006i I2008I) . One re- 
sult is that HCN line emission is enhanced with respect to CO 
in at least some AG N compared with starbursts ( IKohnol l2005t 
IKohno etalll2008b . The high HCN/CO intensity ratios could be 
due to high molecular gas densities near the AGN, and/or high 
HCN/CO abundance ratios that might be due to elevated X-ray 
ionisation and heating rates near the accreting black holes. For 
example, in an early study, ISternberg et al.l (|1994|) concluded 
that the large HCN/CO intensity ratio observed in the nuclear 
(< 100 pc) region in the Seyfert-2 galaxy NGC 1068 indi- 
cates a high HCN/CO a bund ance ratio ~ 10~ 2 near the AGN. 
iLepp & Dalgarnol d 19961) and lUsero etail (120041) argued that el- 
evated HCN/CO ab undance ratios are s i gnatures of high X- 
ray ionisation rates (Malonev et all 19961: Meiierink & Spaansl 



120051: iBoger & Sternberg! 120051: iMeiierink et all 120071) . High 
gas densi ties also likely pla y a role in boosting the HCN line 
intensity (IKrips et al.Ll2008l) . 

In this paper we present high resolution observations of CO 
(2-1) and HCN (1-0) rotational line emissions in the inner re- 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 



gions of NGC 3227, a nearby (D=17Mpc; 1"~ 80pc) Seyfert 
galaxy for which detailed studies have been made of the stellar 
dDavies et alll2006ll2007l) and gaseous dSchinnerer et al.L [20001 : 
iHicks et alll2009l) content of its central regions. We analyze the 
data using large-velocity-gradient (LVG) computations. For this 
galaxy we again find a significant enhancement in the HCN/CO 
intensity ratio close to the active nucleus. The observed nu- 
clear intensity ratio is consistent with optically thick thermal- 
ized emission in dense (> 10 5 cm -3 ) gas. Alternatively, and es- 
pecially in comparison with similar data in NGC 1068 and NGC 
6951, the nuclear emissions could be tracing optically thinner 
emission in which the HCN/CO abundance ratio is large. 



2. Observations 

The analysis in this paper is based on subarcsecond resolution 
observations of the CO(2-l) 230.5 GHz and HCN(l-O) 88.6 GHz 
lines. The CO(2-l) data, for which the beam is 0.6", were previ- 
ously presented bv lSchinnerer et ail d2000l) . New 3 mm HCN(1- 
0) data, with a 0.9"xl.2" beam, were obtained during February 
2009 in the A configuration (760 m baseline) of the six 15- 
meter antennas of the IRAM Plateau de Bure Interferometer. 
The H 12 CN(l-0) line at 88.6 GHz and the H 13 CN(l-0) line at 
86.3 GHz were observed together, using a single polarisation for 
each 1 GHz bandwidth segment. The system temperature was 
80-100 K. Atmospheric conditions were moderate, with winds 
and ~ 5 mm precipitable water vapour. Phase and amplitude 
variations were calibrated out by interleaving reference obser- 
vations of standard calibration sources. The data were processed 
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Fig. 1. Image s of the CO(2-l) line emission (from 
ISchinnerer etall l2000h and HCN(l-O) line emission from 
NGC 3227. Overdrawn is the ellipse tracing the circumnuclear 
ring dDavies et all l2006h . While the CO emission originates 
mostly from the ring, the HCN intensity is far higher in the 
central region. There is an order-of-magnitude change in the 
flux ratio between these regions. 

Table 1. CO and HCN fluxes in NGC3227 



aperture 



Jykms 



Fa 
cai-i 

Jy km s~ 



"//cjvi-o/Fom-i 



1" 0.48 + 0.06 4.4 + 1.3 

2" 1.37 ±0.17 19.5 ± 5.9 

4" 1.86 + 0.29 57.6+ 17.3 

2-4" 0.49 + 0.33 38.5 + 18.3 



114 +UUM 
074 +M33 

-0.023 

033+ 11016 
"■"■"-o.oii 

QQ13+0.024 



w Scaled from the 75mJy given by ISchinnerer et al. d200Ch in a 
8.4"x8.4" aperture and with an uncertainty of 30%. 



and calibrated using the CLIC program in the IRAM GILDAS 
package, and binned spectrally to a resolution of 50 km -1 to 
increase the signal to noise. The kinematics of the HCN(l-0 
line emission in NGC 3227 are analysed by ISani et al ] (120111 
together with several other galaxies. 

It is remarkable that despite the similar beam sizes, the 
CO(2-l) map clearly shows the circumnuclear ring at a radius 
of 140 pc (1.7"), while the HCN map shows emission from 
only the central region, within a few tens of parsecs from the 
nucleus. This can be seen clearly in Fig. Q] and is reflected in 
the line fluxes given in Table. Q] In the nuclear region, which 
we define here to be the central 1" (i.e. r < 40 pc), the flux 
ratio is Fncivi-o/Fco2-i=0.11; on the other hand in the cir- 
cumnuclear ring, for which we take an annulus 2-4", we find 
Fhcni-q[Fco2-i =0.01, an order of magnitude less. Since the 
emission is spatially resolved in both cases, this difference can- 
not be due to differing beam dilution effects. Instead it is likely 
to be due to differences in either molecular abundance - e.g., 
as a result of a distance dependent X-ray ionisation rate - or 
differences in excitation efficiencies - e.g., due to higher nu- 
clear gas densities favoring HCN. (The critical H2 density of 
~ 7 x 10 3 crrT 3 for collisional deexcitation of CO (2-1) is signif- 
icantly lower that the ~ 3 x 10 6 cm 3 for HCN (1-0).) In the next 
section we use large velocity gradient (LVG) models to address 
these issues. 



3. Physical Properties of the Molecular Gas 

3.1. LVG Calculations 

We have constructed a new LVG code and have used it to com- 
pute the HCN (1-0) and CO (2-1) line intensities for a wide 



range of conditions. We calculate the line source functions as- 
suming photo n escape probabil ities from spherical clouds. We 
use the recent lYang et al.l d2010t) data for the excitations and de- 
excitations of the CO rotational level s that are induced by colli- 
sions with H2. For HCN we use the iGreen & Thaddeusl ( fl974h 
collisional data, as updat ed and listed in the RADEX database 
dVan der Tak et al l l2007h . We have cross-checked all of our re- 
sults with the LVG code RADEX, and find excellent agreement: 
in the parameter space we have assessed, the converged line ra- 
tios agree to within 1.5%. We assume that the HCN and CO 
molecules are mixed uniformly, and that the corresponding line 
emissions arise from gas at the same temperature and density. 
The luminosity ratios resulting from the model calculations are 
presented graphically in Fig. [2] covering the parameter space: 
kinetic temperature 30 < T [K] < 300, HCN to CO abun- 
dance ratio 10~ 5 < X HCN /X C0 < 10~ 2 (with X co = 10~ 4 is 
the CO abundance relative to hydrogen), H2 volume gas den- 
sity 10 3 < tin 2 [cm -3 ] < 10 7 , and a ratio of gas-density to 
velocity-gradient, or equivalently column density to linewidth, 
of 5 x 10 17 < N H2 ldV[cm- 2 {kms- 1 )- 1 ] < 5 x 10 23 . 

In each panel in Fig. [2] the parameter-space consists of four 
regimes. The upper right parts correspond to local thermal equi- 
librium (LTE) in the optically thick limit. In this regime the 
transition excitation temperatures T ex of the lines approaches 
the kinetic temperature 7^,,, of the gas. Since the line flux, in 
the Rayleigh- Jeans limit, is F oc J T ex v 2 dv, and an interval dV 
in velocity space is dV = cdv/v, the ratio of two line fluxes 
measured as J F dV (e.g. in units of Jykms -1 as used here) is 

F\ IF2 — (vi / V2) 2 in the optically thick and LTE limit. This cor- 
responds to 0. 15 for the ratio of the HCN(l-O) and CO(2-l) lines 
at 88.6 GHz and 230.5 GHz respectively, which is close to the 
observed ratio of 0. 1 in the nucleus. 

The upper left regions also correspond to LTE because the 
gas densities are high, but in this regime the line optical depths 
are low because the velocity gradients (or line widths) are large. 
For such conditions, the line intensities are linearly proportional 
to the molecular abundances. This behavior is reflected in the 
HCN/CO intensity ratios indicated by the contour values. On the 
left sides of the panels, the intensity ratios decrease linearly with 
the assumed abundance ratio Xhcn/Xco, which ranges from 
10 -2 to 10 -5 in Fig. 

The lower parts of each panel correspond to low densities 
for which the HCN is subthermally excited, leading to relatively 
higher populations in the lowest rotational levels. Again the op- 
tical depth increases from left to right as the column density in- 
creases for fixed line width. 

Some regions of the parameter space may be less physi- 
cally plausible because they assume very large velocity gra- 
dients. Lines of constant velocity gradient are indicated by 
the dashed orange lines in Fig. [2] with increasing gradients 
towards the upper left. For self-gravitating virialised clouds 
AV/R ~ n 1 / 2 where AV is the velocity dispersion, R is 
the cloud radius, and n is the gas density . Treating AV/R 
as a velocity gradient (e.g. iGoldsmitH 12001 1) gives dV/dr ~ 
3.1kms -1 pc -1 y/n H Jl0 4 cm 3 , or typically a few kms 1 pc , 
or a few tens in cases of extreme density. In Fig. |2l the virial 
relation oc (NHjdV) 2 is represented by the dot-dash blue 
line in each panel. Clouds that are unbound or at least partially 
pressure confined could be to the left of this line. 
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Fig. 2. LVG calculations for a 4-dimensional parameter space: the [HCN]/[CO] abundance ratio and gas kinetic temperature are 
given at the top of each panel, while the axes of each panel are the gas volume density and ratio between the column density and 
linewidth (or equivalently between the volume density and velocity gradient). The contours show the expected line ratio based on 
the emitted fluxes in Jykms 1 . The locus of representing possible parameters for the nucleus and ring of NGC 3227 are drawn in 
red and green respectively. The dashed orange lines are curves of equal velocity gradient (in km s pc~' ), indicating which regions 
of the parameter space are physically plausible. In particular, the locus for self-gravitating virialised clouds is represented by the 
dot-dash blue line. 



3.2. Analysis for NGC 3227 

Curves representing the possible parameter ranges for the mea- 
sured flux ratios ^hcnx-qI^coi-x for NGC 3227 have been over- 
drawn for the nucleus (red lines) and circumnuclear ring (green 
lines), as defined in Sec. [2] If we assume that the clouds are 
self-gravitating, we are restricted to the points where the red and 
green lines intersect the dashed blue line. These show that if the 
clouds in the ring have a similar density to the nucleus, then 
the HCN abundance must be about 2 orders of magnitude lower. 
Alternatively, the density may differ by up to 2 orders of mag- 
nitude if the abundances are similar. We discuss this further in 
Sec.lBH 

In the nucleus, the flux ratio Fhcni-o/Pco2-\ =0. 1 1 is remark- 
ably close to the LTE optically thick limit, and its locus includes 
a contour around this region. However, depending on the HCN 
abundance, the gas density at which this occurs can vary from 
Hh 2 ~ 10 4 cm~ 3 at the highest abundance to ~ 3 x 10 5 crrT 3 at 
the lowest abundance we have considered, with column density 



Nh^ /dV > 10 22 cirT 2 /(km s ). By considering the mean volume 
density of the HCN emitting region, and putting a limit on a re- 
alistic filling factor, we can restrict this range further. 

The first step is to estimate the volume of the emitting re- 
gion. This can be done because the HCN emission is marginally 
resolved. A detailed estimate of the intrinsic size - via dynami- 
cal model ling, and taking in to account emission from the ring - 
is given in Sani et all d201 lh . These authors show that the diame- 
ter is 0.54", corresponding to 45 pc, and the scale height is 6 pc. 
Taking this as an indication of the size along the line of sight sets 
the volume of the emitting region. 

The second step is to estimate the mass. We have done this 
in several ways since they are all uncertain. 

1. The LVG calculation yields a mass directly under the as- 
sumption that the ratio of the observed linewidth to the 
linewidth of an individual cloud is tracing the number of 
clouds, such that N, ot = N c i m ,d(6v b s /dv c i ou d). For condi- 
tions corresponding to XhcnIXqo = 10~ 2 , T= 300 K, 
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n Hl = 10 5 - 5 crrT 3 and N H JdV = 10 22 cm- 2 (kms" 1 )- 1 (see 
Sec. 13.3b . the observed HCN(l-O) flux leads to a mass of 
1.8 x 1O 6 M . 

2. The standard method to estimate the mass is from the CO lu- 
minosity. We have done this directly from the CO(2-l) line 
flux in Table Q] using a conversi on factor a - 4.3 wh ich in- 
cludes a correction for helium dTacconi et all 120081) . This 
yields 3.3 x 1O 6 M . 

3. A similar conversion f o r the HCN line that has been cal- 
ibrated by iKrips et alJ d2008l) for AGN is M Hl /L HCN ~ 
10 M (Kkm s" 1 pc 2 )" 1 . This yields a mass of 6 x"l0 6 M . 

4. As a final check, we use the dynamical mass derived from 
the HCN kine matics. By fitting models to account for the 
beam smearing. lSani et al.l (1201 ll) find M dyn = 5.6 x 10 7 M . 
For a nom inal 10% gas fract ion expected in local disks and 
starbursts dHicks et all 12.009b . this would yield a gas mass of 
6x 1O 6 M . 

These estimates are all the same order of magnitude, and 
suggest that the gas mass in the central arcsec is of order 
4 x 1O 6 M . Hence we can estimate the mean density to be 
^ 6 x 10 3 cirT 3 . Comparing this to the cloud densities 
above yields volume filling factors in the range 1-0.01. In this 
range, a lower filling factor is more physically plausible, which 
would tend to favour the solutions with higher cloud densities. 
Fig.|2]shows these have either higher temperature or less extreme 
HCN abundance. 

3.3. Comparison to NGC 6951 and NGC 1068 

NGC 1068 and NGC 6951 are two other galaxies for which 
the HCN(1-0)/CO(2-1) ratio has been measured on comparable 
~ 100 pc scales. We use flux densities reported by IKrips et al] 
(120071) for the nuclear region (denoted 'C in their Table 1) of 
NGC 6951; and also the values for the circumnuclear disk of 
NGC 1068, as the sum of the red and blue channels reported 
in Table 3 of lUsero et al.l (120041) . These yield line ratios (for line 
fluxes in Jy km s" 1 ) of 0.37+0.05 and 0.2 14 + 0.002 respectively, 
and are denoted by the solid magenta lines on Fig. [2] These 
lines appear almost exclusively in the panels corresponding to 
the highest HCN abundance we have considered, Xhcn/Xco = 

io- 2 . 

In contrast to NGC 3227, in which the line emission appears 
to be optically thick, the loci of the magenta lines for NGC 1068 
and NGC 695 1 are toward the optically thin (left) side of the pan- 
els. Despite this, it is notable that there are regions of the param- 
eter space where the contours corresponding to all 3 objects lie 
close together, running from lower left to upper right. The region 
extends from «h, = 10 4 crrT 3 and Nh 2 /dV = 10 19 cm _2 /(kms _1 ) 
to nn 2 = 10 6 cm 3 . It is precisely because one can attribute 
the observed line ratios - with different optical depths for the 
3 galaxies - to similar physical properties of the gas in all these 
3 objects that this locus is appealing. 

Why this occurs can be seen in Fig. |3]which shows the op- 
tical depths r for the HCN(l-O) and CO(2-l) transitions. The 
gas properties of both these panels correspond to the bottom left 
panel in Fig. |2](300K and Xhcn/Xco = 10~ 2 ), and cover the 
same range of density and velocity gradient. These plots show 
clearly the characterisation of the different regions: in the lower 
half the HCN(l-O) line is optically thick because the density is 
low enough that it is sub-thermal; above the critical density, the 
line is in LTE and thus optically thin at low columns and op- 
tically thick at high columns. The locus where all the contours 
for the 3 galaxies are close together and parallel follows approx- 



Xhcn/Xco= 1 °" 2 . X co =10" 4 , T=300K 




18 19 20 21 22 23 18 19 20 21 22 23 

log N H2 /dV (cm -2 / km s" 1 ) log N H2 /dV (cm -2 / km s"') 



Fig. 3. Optical depth t of the HCN(l-O) line (left) and CO(2- 
1) line (right) for gas properties corresponding to the bottom left 
panel in Fig. [2] Darker regions correspond to lower optical depth. 
Overplotted are the contours for the 3 galaxies corresponding to 
their HCN(1-0)/CO(2-1) ratio. 



imately the boundary where the HCN(l-O) line becomes opti- 
cally thick. Here, a small change in physical conditions (column 
or density) can result in the HCN(l-O) emission switching from 
optically thin to optically thick. 

This regime is, however, also associated with very large ve- 
locity gradients. It is dV/dr ~ 10 4 kms _1 pc _I at T = 30 K, 
but reduces as the tempera ture increases. Velocity gr adients were 
not dis cussed explicitly bv lSternberg et alJ ( 1 19941) or lUsero et al.l 
(12004 in their T = 50 K LVG calculations for NGC 1068. But 
their analyses also associate the observed properties with sim- 
ilarly extre me velocity gradients. Indeed, one of the main con- 
clusions of lSternberg et al.l (11994 was that X HC n/X C o S 10" 2 
in NGC 1068. For the temperature they considered, this would 
lead to dV/dr ~ 10 4 kms -1 pc -1 (matching the top left panel of 
Fig. |2] here). However, our LVG calculations shows that dV/dr 
is reduced as both the temperature and density increase. When 
considering all 3 galaxies together, the smallest - and therefore 
arguably the most physically plausible - value in the parameter 
space we have covered is dV/ dr ~ 1 00 km s~ 1 pc~ 1 at T = 300 K 
and «h, ~ IO 5 ' 5 cm 3 . This location is not far from the bound- 
ary of the optically thick LTE regime discussed previously, but 
due to the high velocity gradient represents clouds that are either 
pressure confined or unbound. 

Interestingly, there is evidence in NGC 1068 from recent 
Herschel observations with PACS of high rotational CO tran- 
sitions, for a significant mass of molecular gas in the central 
~ lOOpc at te mperatures of 100K and 400 K and densities of 
~ 10 6 5 cm 3 dHailev-Dunsheath et ail 1201 ll) . Similarly, in an 
analysis of various HCN, HCO + and CO isotope t ransitions in 
the central lOOpc of NGC 1068. lKrips et al.l d201 ll) also argued 
in favour of warm (T> 200 K) gas. However, they also concluded 
that the density is of order nn 2 ~ 10 4 cm 3 . Our LVG calcula- 
tions indicate that such densities are associated with very high 
velocity gradients for the observed line ratio, which we consider 
physically unlikely. This, combined with a comparison of the 
cloud density to mean density estimated in Sec. l3.2l has led us to 
favour the higher density solution with more moderate velocity 
gradient. 

3.4. X-ray lonisation rate in NGC 3227 

Our LVG calculations presented above suggest that as one 
alternative the high HCN(1-0)/CO(2-1) ratio in the central 
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~100pc of NGC 3227 could be due to an exceptionally 
high HCN abundance, an interpretation supported by the 
even higher HCN/CO intensity ratios observed in the nu- 
clei of NGC 1068 and NGC 6951. It is possible that the 
high HCN abundances are associated with elevated X-ray 
ionisation and/or heating rates near the AGN. Theoretical 
investigations of X-ray (or cosmic -ray) driven chemistry 
show that the (steady-state) molecular abundances depend 



isation rate f ([Rrolik & Kallma 


d 1983: Leoo & Dalprno 


19961 Malonev et all 1 19961: 


Meiierink & Spaans, 2005; 


Boger & Sternberg, 2005; Meiierink et alll2007l). Here we adopt 



the notation of iBoeer & Sternberg! (120051) . normalising ( to 
10~ 17 s _1 to give the parameter ra/f-17. The density n — nn + «h 2 
refers to the total atomic plus molecular hydrogen density. In 
the following analysis, we first calculate the ratio n/g-n in the 
nucleus of NGC 3227 and then compare it to ratios predicted by 
models of X-ray irradiated gas. 



3.4.1. n/£-n in NGC 3227 

To estimate the X-ray ionisation rate in the central ~100pc 
of NGC 3227 we need to know the intrinsic spectral energy 
distribution from the AGN. We adopt the SXPL model of 
iMarkowitz et alj d2009l) in which both the hard and soft com- 
ponents of the X-ray flux are modelled with power laws: 

7Y ph = 0.0040(£/A:eVr 3 - 35 + 0.0067(£//teVr L57 

where E is the photon energy in keV and N p h is the photon flux 
in units of phkeV 1 cirT 2 s _1 . Making the usual assumption that 
the primary ionisation rate of hydrogen is negligible compared 
to the secondary ionisation rate, we t hen calculate the resu lting 
ionisation rate f using equation A4 of Mal onev et al.l d 1996b : 



c 



ctUE) F(E)dE. 



Here N sec = 28 (as given by iMalonev et al.llT996b is the num- 
ber of secondary ionisations per keV of primary photoelectron 
energy assuming a mean-energy per ion-pair of 37.1 eV for en- 
ergy deposition in a molecular hydrogen gas (Dal garno et al.l 
Il999h : cr pa (£') is the absorption cross section per H nucleus, 
for which we adopt th e broken power-law fit in equation A5 of 
IMalonev et al.l (Il996l) : and F(E) is the incident flux in units of 
ph keV~ 1 . The photoionisation is dominated by photons with en- 
ergies for which t ~ 1 . We therefore take the limits of the in- 
tegral to be E max = lOOkeV and E m i n as the energy at which 
the optical depth due to photoelectric absorption is t(E) = 
1, ignoring attenuation above this limit. T he energ y at w hich 
t( E) = 1 is interpolated from Tabl e 9.3 of lSewardi (l2000l) for 
the iMorrison & McCammonl (Il983l) model, for a given column 
density. And the column density is assumed to be proportional 
to distance from t he AGN up to a maximum of 3 x 10 23 cm -2 
dHicks et al.ll2009h at 30 pc. 

We have evaluated the integral at two distances: 18pc, an 
area weighted mean distance from the AGN to the HCN emit- 
ting gas that corresponds to the nuclear region; and 140 pc cor- 
responding to the distance of the circumnuclear ring. We find 
£ 18pc = 3.6xl0 _13 s _1 andfi4o pc = 3.5 x 10~ 15 s" 1 , about a factor 
100 less primarily due to the distance related geometrical dilu- 
tion of the incident X-ray flux. The inferred ionisation rates are 
much larger than the typical ionisation rates in Galactic clouds. 

Adopting a characteristic density «h, ~ 10 5 ' 5 crrT 3 from our 
LVG analysis in Sec. 13.31 yields n/£-n ~ 10 in the nuclear re- 
gion of NGC 3227. We note that even if the gas density were 
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Fig. 4. Abundance ratio [HCN]/[CO] as a function of ra/f-n, the 
ratio of the local gas volume density to ionisation rate, for 100 K 
gas, showing the transition from the high- to low-ionisation 
phase. 



an order of magnitude higher (leading to proportionally higher 
n/^-n), the physical conditions will be well within the high- 
ionisation phase for the molecular chemistry. 

3.4.2. Models of abundance ratio as a function of «/f_n 

The model computations of iBoger & Sternberg! d2005l 120061) 
show that gas can exist in a high (low) ionisation phase for small 
(large) values of n/g-n, with a density ratio of tt/f-n ~ 10 3 
marking the cross-over between the two regimes. The HCN/CO 
abundance ratio can become large > 10~ 3 in the high-ionisation 
phase even at low gas temperatures, due to the large densities 
of atomic and ionic carbon. Here we have re-run these calcula- 
tions using the same elemental gas-phase abundances, and with 
a slightly updated reaction set. Fig. |4] shows the resulting abun- 
dance ratios as a function of for 100 K gas. At very low 
n/^-n < 100 corresponding to that in the nuclear region of 
NGC 3227 the HCN/CO abundance ratio approaches ~ 10 3 . 
For comparison, iLepp & Dalgarnol (1 19961) find a peak value of 
5 x 10~ 4 in their calculation. These models track the abundances 
of numerous molecules across many orders of magnitude and 
show that at low n/f-n, the HCN/CO abundance ratio is raised 
several orders of magnitude above that typically expected, to a 
level at which it approaches - and, given the uncertainties in 
such models, is commensurate with - that implied by the LVG 
mo dels. 

lHarada et al.l (1201 Oh have shown that at elevated tempera- 
tures > 300 K, rapid hydrogenation of CN to HCN can increase 
the HCN/CO abundances further still, and they comment that 
such warm gas may be present in the X-ray heated gas near 
AGN. We caution however, that their higher temperature mod- 
els for which HCN/CO is largest, correspond to n/f-17 = 10 4 ' 5 
which is significantly larger than the value we are invoking as 
characteristic for the nucleus of NGC 3227. 

The models show that a combination of X-ray ionisation and 
heating do yield high HCN abundances, although not yet quite 
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as high as HCN/CO~ 10~ 2 as implied by the LVG analysis in 
Sec. l3.3l for the nucleus. Further chemical modeling is required, 
but a high HCN/CO intensity ratio in the nucleus due to ele- 
vated X-ray ionisation rates appears plausible. In this picture, 
the lower HCN/CO intensity ratio in the ring may simply reflect 
the lower ionisation rate there, and not just a lower gas density 
in this circumgalactic environment. 

4. Summary and Conclusions 

We present an LVG analysis of high-resolution observations of 
CO (2-1) and HCN (1-0) line emissions in the central regions of 
the Seyfert galaxy NGC 3227. We find that 

- The HCN(1-0)/CO(2-1) ratio in NGC 3227 is an order of 
magnitude higher in the central 80 pc than in the circumnu- 
clear ring at a radius of 140pc. NGC 6951 and NGC 1068 
have similarly high published ratios in their central ~ 100 pc. 

- The nuclear HCN/CO intensity ratio in NGC 3227 may re- 
flect optically thick line emission in dense gas with only a 
weak constraint on the HCN/CO abundance ratio. However, 
our LVG calculations also indicate that the high nuclear ra- 
tios in all three of these galaxies are more consistent with 
a single set of physical properties corresponding to warm 
~ 300 K, dense 10 55 crrT 3 gas, in which the emission lines 
are optically thinner, but in which the HCN/CO abundance 
ratio is very large ~ 10~ 2 . For these conditions the velocity 
gradients are dV/dr ~ lOOkms -1 , but would increase sig- 
nificantly at lower temperatures or densities. Most likely the 
clouds are gravitationally unbound. 

- The X-ray ionisation rate at radii less than ~ 20 pc may ex- 
ceed 10~ s , and could plausibly lead to high HCN abun- 
dances in molecular gas in the high-ionisation phase where 
the ratio of the gas density to the X-ray ionisation rate is 
small. 
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